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Abstract 

In order to understand whether, and to what extent, spectral 
representation can effectively highlight the nonlinear interac- 
tion among different scales, it is necessary to consider the state 
that precedes the onset of instabilities and turbulence in flows. 
In this condition, a system is still stable, but is however sub- 
ject to a swarming of arbitrary three-dimensional small per- 
turbations. These can arrive any instant, and then undergo a 
transient evolution which is ruled out by the initial- value prob- 
lem associated to the Navier-Stokes linearized formulation. The 
set of three-dimensional small perturbations constitutes a sys- 
tem of multiple spatial and temporal scales which are subject 
to all the processes included in the perturbative Navier-Stokes 
equations: linearized convective transport, linearized vortical 
stretching and tilting, and the molecular diffusion. Leaving 
aside nonlinear interaction among the different scales, these fea- 
tures are tantamount to the features of the turbulent state. 

We determine the exponent of the inertial range of arbitrary lon- 
gitudinal and transversal perturbations acting on a typical shear 
flow, i.e. the bluff-body wake. Then, we compare the present 
results with the exponent of the corresponding developed turbu- 
lent state (notoriously equal to —5/3). For longitudinal pertur- 
bations - i.e. perturbations in the plane of the basic flow which 
is two-dimensional - we observe a decay rate of —3 in the in- 
ertial range, typically met in two-dimensional turbulence. For 
purely three-dimensional perturbations, instead, the energy de- 
creases with a factor of —5/3. If we consider a combination 
of longitudinal and transversal perturbative waves, the energy 
spectrum seems to have a decay of —3 for larger wavenumbers 
(k £ [50, 100]), while for smaller wavenumbers (k e [3, 50]) the 
decay is of the order —5/3. We can conclude that the value 
of the exponent of the inertial range has a much higher level of 
universality, which is not necessarily associated to the nonlinear 
interaction. 

Introduction 

One very popular notion in the phenomenology of turbulence 
(in the sense of Kolmogorov 1941) is that a power-law scal- 
ing with an exponent close to —5/3 is observed for the en- 
ergy spectrum over a very substantial range of a few decades 
of wavenumber, this range being called the inertial range. For 
extensive collections of laboratory and numerical experimental 
results, see for instance (TJ [2j . In fact, it is a common crite- 
rion for the successful production of a fully developed turbulent 
field, either in the laboratory or in numerical simulations, to 
verify that the power spectrum has such a scaling in the inertial 
range. 

The set of arbitrary three-dimensional small perturbations con- 
stitutes a system of multiple spatial and temporal scales which 
are subject to all the processes included in the Navier-Stokes 
equations, leaving aside the nonlinear interaction among the dif- 



ferent scales. If it were possible to observe such a system in a 
temporal window and obtain the instantaneous power spectrum, 
it would be possible, among others, to determine the exponent 
of the inertial range of the arbitrary perturbation, and to com- 
pare it with the exponent of the corresponding developed tur- 
bulent state. Two possible situations can therefore appear: (a) 
the exponent difference is large and, as such, it is a quantita- 
tive measure of the nonlinear interaction in spectral terms; (b) 
the difference is small. This would be even more interesting, 
because it would indicate a higher level of universality on the 
value of the exponent of the inertial range, not necessarily asso- 
ciated to the nonlinear interaction. 

We propose building temporal observation window for the tran- 
sient evolution of a large number (order of 10 2 ~ 10 3 ) of arbi- 
trary small 3D (oblique) perturbations acting on a typical shear 
flow. In particular, we can take advantage of a recently - nu- 
merically obtained - set of solutions yielded by the initial- value 
problem applied to 3D perturbations of a plane bluff-body wake 
(3JH). These solutions have revealed the existence of many 
different kinds of transient behaviour, not all of which is triv- 
ial. If these transients, obtained in association with arbitrary 
initial conditions, are injected in a statistical way into the tem- 
poral observation window, we can obtain a close representation 
of the perturbation state that precedes the onset of instability- 
turbulence. In this preliminary work, we consider an ensemble 
of transients of asymptotically stable waves - this is because 
we want to avoid any temporal divergence - and build the en- 
ergy spectrum by freezing each wave at the moment it reaches 
a constant value of the temporal growth rate. Here, the growth 
is always negative since all waves are asymptotically damped 
in time. When this happens, the wave may be considered out 
of its transient and in the asymptotic condition. When build- 
ing the power spectrum, the energy of each wave is normalized 
over its initial energy. In this way, we are considering the spec- 
tral dynamics of a sort of white noise perturbation made up of 
asymptotically stable waves. 

Formulation 

The energy spectrum behaviour is studied using the initial- value 
problem formulation. The base flow is approximated at a fixed 
longitudinal station, xq = 10, through an analytical expansion 
solution | 5 | of the Navier-Stokes equations (see Fig. [TJ. The 
Reynolds number is set to a value of 40, in order to consider sta- 
ble evolutive configurations. The viscous perturbative equations 
are written in terms of the vorticity and the transversal velocity 
1 6 1 and then transformed through a Laplace-Fourier decompo- 
sition (3][4) in the plane (x,z) which is normal to the base flow 
plane (x,y), 
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Figure 1: Perturbation geometry scheme, mean flow U in the 
wake at Re = 40 and xo = 10 (blue curve), symmetric and asym- 
metric initial conditions in terms of v(f = 0, y) (red curves). 
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The transversal velocity and vorticity components are indicated 
as v and u) v respectively, while f is defined through the kine- 
matic relation T = d x CQ z — 3 ; tt>i that in the physical plane links 
together the perturbation vorticity components in the x and z di- 
rections (Cflf and co ; ) and the perturbed velocity field. Equations 
{2} and ^ are the Orr-Sommerfeld and Squire equations re- 
spectively, from the classical linear stability analysis for three- 
dimensional disturbances in the phase space. We define k as 
the polar wavenumber, a,- = kcos(§) as the wavenumber in x 
direction, y = ksin(§) as the wavenumber in z direction, as 
the angle of obliquity with respect to the physical plane, and oc, 
as the spatial damping rate in x direction (see the perturbation 
scheme in Fig. [TJ. The measure of the perturbation growth can 
be defined through the disturbance kinetic energy density in the 
plane (a,y): 
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where u and w are the streamwise and spanwise components of 
the perturbation velocity, respectively, while ly^ is the exten- 
sion of the spatial numerical domain. The amplification factor 
G(t) can be introduced in terms of the normalized energy den- 
sity 



G(t;a,y) = 
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Figure 2: (a) The amplification factor G, asymmetric initial 
condition, intermediate (xo = 10, solid curves) and far field 
(xo = 50, dashed curves) wake configurations. The periods 
tinier, 1 f ar are the periods of the modulation visible on G, in 
the intermediate and far field, respectively, (b) The amplifica- 
tion factor G and (c) the temporal growth rate r as function of 
time, symmetric initial condition, k = 0.5, 1, 1.5,2. 
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is introduced in order to evaluate both the early transient as well 
as the asymptotic behaviour of the perturbations. 

Examples, in terms of G and r, of significant transient behaviour 
and asymptotic fate of the three-dimensional perturbations 1 3 1 
are shown in Fig. [2] In part (a) the amplification factor G 
of an asymmetric wave is shown for two typical intermediate 
(xo = 10) and far (xo = 50) wake configurations. For xq = 10 
a local maximum, followed by a minimum, is visible in the en- 
ergy density, then the perturbation is slowly amplifying and the 
transient can be considered extinguished only after hundreds of 



time scales. For xo = 50 these features are less marked, but still 
present. For both configurations, the function G shows a mod- 
ulation in time in the first part of the transient (%f ar =100 and 
timer = 35), which is always observed in the case of asymmet- 
ric longitudinal or oblique instability waves. In parts (b) and (c) 
of Fig. 12 the amplification factor G and the temporal growth 
rate r of purely three-dimensional symmetric perturbations are 
reported, respectively. Such waves may become damped by in- 
creasing their wavenumber (k = y). Before the asymptotic sta- 
ble states are reached, these configurations yield maxima of the 
energy density (e.g. when /: = 1 .5, G ~ 3 at ? ~ 30) in the tran- 
sients. This trend is also typical of oblique and longitudinal 
waves, and it can be considered a general feature in the context 
of the stability of near parallel shear flows. 

Results 

Computations to evaluate the energy spectrum are made, at a 
fixed wavenumber, by integrating the equations forward in time 
until the temporal growth rate r asymptotes to a constant value, 
i.e. when the transient can be considered as extinguished. We 
consider a range of values for the polar wavenumber k in the in- 
terval [0.05, 100], and three angle of obliquity (|> = 0,Jt/4,Jt/2. 
For simplicity's sake, the spatial damping rate oc, is here taken 
equal to 0, thus no spatial damping is considered in the lon- 
gitudinal direction. We account for symmetric and asymmet- 
ric initial conditions (see red curves in Fig. [T} in terms of the 
transversal velocity v, while the transversal vorticity (h y is ini- 
tially equal to zero. 

The normalized energy density G at the asymptotic state (blue 
and black circles for symmetric and asymmetric initial con- 
ditions, respectively) and the initial energy density e(t = 0) 
(blue and black triangles for symmetric and asymmetric ini- 
tial conditions, respectively) are shown - as function of the 
polar wavenumber k - in parts (a) and (b) of Fig. [3] for a 
longitudinal ((]) = 0) and a transversal (<]) = jt/2) perturbative 
wave, respectively. In Fig. [5J; we sum, at a fixed wavenumber, 
asymmetric and symmetric perturbations with angle of obliq- 
uity (]) = 0,7r/4,7r/2, and then we calculate the energy density 
of the resulting perturbation obtained by combining six differ- 
ent waves. We report —5/3 and —3 slopes with red and black 
curves, respectively. 

In the case of longitudinal waves (two-dimensional base and 
perturbed flow configurations, see Fig. [3ji), the normalized en- 
ergy density G at asymptotic state has a decay of —3 in the iner- 
tial range (k £ [2, 100]) for both symmetric and asymmetric per- 
turbations (see blue and black circles). For smaller wavenum- 
bers, instead, the decay is faster. For purely transversal waves 
(Fig. [3p), we observe a decay of —5/3 in the inertial range 
(k £ [2, 100]) for both symmetric and asymmetric perturbations 
(blue and black circles). Longer waves have a deeper decay, 
reaching a maximum of energy at about k = 0.5. If we sum sym- 
metric and asymmetric waves with different angles of obliquity 
((|> = 0, Jt/4, Jt/2), the normalized energy density G (blue cir- 
cles) of the resulting perturbation at the asymptotic state has a 
decay of —3 for k £ [25, 100], while the decay is of the order 
— 5/3 for k £ [2,25]. For longer waves, the decay is faster and a 
maximum of energy can be found for k £ [0.1,0.5]. 

From Eq. l|4j and recalling that the transversal vorticity d) v is 
initially imposed equal to zero, one can observe that the initial 
energy density has a well-defined decay, given by: 
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Figure 3: Spectrum of the amplification factor of a collec- 
tion of stable perturbation waves, Re = 40. The normalized 
energy density G at the asymptotic state (circles) and the initial 
energy density eit = 0) (triangles). The energy of each wave 
in asymptotic conditions is normalized over the value owned 
at the initial instant. Blue and black symbols for symmetric 
and asymmetric initial conditions, respectively, (a) (j) = 0, (b) 
<|) = 3l/2, and (c) sum of symmetric and asymmetric perturba- 
tions with angle of obliquity <|) = 0,7l/4,7t/2. Red and black 
curves indicate —5/3 and —3 slopes, respectively. Note that the 
intermediate waves are stable and highly damped perturbations, 
that reach the asymptotic state after a long lapse of time where 
they progressively lose their kinetic energy. This explain the 
sharp fall of energy density at the transition between the long 
and the intermediate wave ranges. 
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This analytical behaviour is confirmed for the three cases we 
considered in Fig. [3] In parts (a) and (b) blue and black triangles 
(for symmetric and asymmetric initial conditions, respectively) 
represent the initial energy density and show a constant value 
for large wavenumbers (k € [2, 100]), while a decay of about 
—2 is visible for longer waves (k e [0.05, 2]). It should be noted 
that eit = 0) does not depend on the perturbation obliquity, thus 
the energy density values at t = are exactly the same for both 
(|) = and (]) = 7T./2. In case (c), where six different waves are 
summed and then the energy density of the resulting perturba- 
tion is computed, e(t = 0) (blue triangles) still displays a trend 
in agreement with Eq. 0. 
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Conclusions 

This is a preliminary study of the behaviour of a collection of 
transient waves seen throughout their power spectrum. To avoid 
the complexity introduced by the temporal divergence associ- 
ated to asymptotically unstable waves, we considered an en- 
semble of transients of asymptotically stable waves. Then, we 
have built the energy spectrum by freezing each wave at the 
moment it reaches a constant (negative) value of the temporal 
growth rate. When this happens, the waves are dynamically out 
of their transient and in the asymptotic condition. In the power 
spectrum, the energy of each wave is normalized over its own 
initial value. In practice, it is the same as considering the spec- 
tral dynamics of a white noise perturbation, a kind of model for 
the swarm of small perturbations that affects any system in a 
linear way. 

We observe that, whether the waves are aligned with the base 
sheared flow or not, the energy of the intermediate range of 
wavenumbers in the spectrum decays with the same exponent 
(—5/3 for 3D oblique waves, —3 for 2D aligned waves) that 
is observed in the spectrum of the velocity fluctuation of fully 
developed turbulent flows. Where the nonlinear interaction is 
considered dominant. 

At the moment, we can conclude by observing that the spectral 
power-law scaling of intermediate/inertial waves (with an expo- 
nent close to —5/3 in 3D and to —3 in 2D) is a general dynam- 
ical property which encompasses the nonlinear interaction. In 
other words, it seems to us that the strength of the nonlinearity 
has been overestimated in the determination of a few turbulence 
properties. 
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